giimq 


Eleventh ARM Science Team Meeting Proceedings, Atlanta, Georgia , March 19-23, 2001 


Cloud Overlapping Detection Algorithm Using Solar and IR 
Wavelengths with GOES Data Over ARM/SGP Site 


K. Kawamoto 

Virginia Polytechnic Institute and State University 
Blacksburg , Virginia 

P. Minnis and W. L. Smith , Jr. 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 


Introduction 

One of the most perplexing problems in satellite cloud remote sensing is the overlapping of cloud layers. 
Although most techniques assume a one -layer cloud system in a given retrieval of cloud properties, 
many observations are affected by radiation from more than one cloud layer. As such, cloud overlap can 
cause errors in the retrieval of many properties including cloud height, optical depth, phase, and particle 
size. A variety of methods have been developed to identify overlapped clouds in a given satellite imager 
pixel. Baum et al. (1995) used CO 2 slicing and a spatial coherence method to demonstrate a possible 
analysis method for nighttime detection of multi-layered clouds. Jin and Rossow (1997) also used a 
multi- spectral CO 2 slicing technique for a global analysis of overlapped cloud amount. Lin et al. (1 998) 
used a combination infrared (IR), visible (VIS), and microwave data to detect overlapped clouds over 
water. Recently, Baum and Spinhime (2000) proposed a 1 .6 and 1 1 pm bispectral threshold method. 
While all of these methods have made progress in solving this stubborn problem, none have yet proven 
satisfactory for continuous and consistent monitoring of multi-layer cloud systems. It is clear that 
detection of overlapping clouds from passive instruments such as satellite radiometers is in an immature 
stage of development and requires additional research. Overlapped cloud systems also affect the 
retrievals of cloud properties over the Atmospheric Radiation Measurement (ARM) domains (e.g., 
Minnis et al. 1998) and hence should be identified as accurately as possible. To reach this goal, it is 
necessary to determine which information can be exploited for detecting multi-layered clouds from 
operational meteorological satellite data used by ARM. This paper examines the potential information 
available in spectral data available on the Geostationary Operational Environmental Satellite (GOES) 
imager and the National Oceanic Atmospheric Administration (NOAA) Advanced Very High 
Resolution Radiometer (AVHRR) used over the ARM Program’s Southern Great Plains (SGP), and 
North Slope of Alaska (NSA) sites to study the capability of detecting overlapping clouds. 

Data 


This study uses daytime half-hourly GOES-8 4-km data from channels 1, 2, 4, and 5, at VIS 0.65 pm, 
solar- IR 3.9 pm, IR 10.8 pm, and split-window (WS) 12.0 pm wavelengths, respectively, and the cloud 
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properties over the ARM SGP central facility (CF) derived from these radiances by Minnis et al. (2001) 
for all of 1998. The ARM 35-GHz radar data taken over the CF are used to provide cloud boundary data 
to determine the presence of single- and multi-layered clouds during each 10 - minute period centered on 
a given GOES image time. The average GOES radiances and cloud properties for a 0.3° box centered 
on the CF are compared to the cloud boundary data from the corresponding radar data. 

Methods and Results 

Two methods for detecting multi-layered clouds are proposed. The first one relies on the brightness 
temperature difference between channels 4 and 5 (hereafter BTD 45 ). The basic assumption for this 
approach is that for large cloud optical depth x, BTD 45 should be small, nearly equal to zero, if the cloud 
is a single layer because the emissivities in both channels approach unity. On the other hand, if BTD 45 
is large and x is large, a thin cloud probably overlaps a lower cloud. Thus, BTD 45 should be an indicator 
of overlapping clouds when the x is large. To illustrate this situation. Figure 1 shows the results of a 
numerical simulation using the parameterizati ons of Minnis et al. (1998). The values of BTD 45 are 
plotted as functions of total optical depth (OD) and effective cloud particle size for single and 2- layer 
cases. The single cloud is ice layer with OD between 0.5 and 9.5 with an effective diameter De = 1 8, 30, 
68 , and 123 pm The two-layer system consists of a low cloud x = 10 and an effective droplet radius 
r e = 10 pm overlaid by the high cloud used in the single-layer case. Cloud-top temperatures of 245 and 
275 K are used for the high and low clouds, respectively. The figure demonstrates that the values of 
BTD 45 for the two-layer cloud system are greater than their single- layered counterparts for a given value 
of total OD. In this case, values of BTD 45 exceeding about 0.6K for x > 8 should indicate the presence 
of multi-layered clouds. Although a great variety of cloud temperature, OD, and particle size 
combinations will occur, BTD 45 values should follow a pattern similar to that in Figures 1 , 2, and 3 plot 
the observed BTD 45 values as a function of the retrieved OD for both single- and multi-layered clouds 
having BTD 45 > 0.5 K and x > 10 correspond to multi-layered clouds. For smaller BTDs, only 
21 percent of the clouds are multi-layered. If a BTD 45 of 1 K is used as the threshold, 85 percent of the 
clouds are multi-layered for larger BTD 45 , but 28 percent of the multi-layered clouds would be identified 
as sing le-layered. At lower ODs, BTD 45 is generally greater than that from the single- layer clouds, but 
no clear threshold values are evident for selecting single- versus multi-layered clouds. 

For clouds identified as ice (Figure 3), only 60 percent of the pixels are clearly multi-layered using a 
threshold of BTD 45 = 0.5 K for the optically thick clouds. Many of the optically thicker clouds have 
BTD 45 > 1 .0 K. If the threshold of 1 .0 K is used for the ice clouds, 82 percent of the clouds with larger 
BTD 45 are multi-layered and 36 percent of the clouds with smaller BTD 45 are single layered. The 
differences between Figures 2 and 3 may be due to the low density of ice clouds and their IR scattering 
properties. Because the phase discrimination is based on models of cloud particle sizes, the ice cloud 
emissions probably dominate the clouds identified as ice, while those determined to be water most likely 
have only a thin layer of ice cloud over the low-level water cloud. The water clouds are generally dense 
so the radiating temperatures for both channels in the single-layer case are probably very close. 

Optically thin ice clouds can be several kilometers thick with vertically dependent microphysical 
properties. Optically thicker ice clouds also may be very thick physically with a layer of low- density 
cirrus at the top of the cloud. These low- density clouds or cloud layers could easily produce BTD45 = 

1 K, even for optically thick clouds because of their vertical structure. Thus, detecting multi- layering 
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Figure 3. Same as Figure 2, except for clouds identified as ice phase only. 

for clouds identified as ice is more difficult than for water clouds using this approach. The second 
method uses particle size for detecting overlapping clouds. The method applied by Minnis et al. (2001) 
assumes that the observed cloud in a given imager pixel is a single-layered cloud and uses the 3.7- and 
1 1 -pm brightness temperature difference between channels 3 and 4 (hereafter BTD 34 ) to determine 
phase and effective particle size. For a given value of x, BTD 34 is smallest for large ice crystals and 
greatest for small droplets. From the smallest water droplet radius, droplet size generally increases as 
BTD 34 decreases. Conversely, starting with the largest ice crystal, effective diameter (D e ) decreases as 
BTD 34 increases. The value of BTD 34 may correspond to both large water droplets and small ice 
crystals for intermediate cases. Thus, if a relatively thin ice cloud overlaps a water cloud, the retrieved 
value of effective radius (r e ) or D e may be either very large or extremely small, respectively, depending 
on the OD of the upper level cloud. Motivated by this idea, effective particle sizes were plotted against 
the retrieved OD values for single and multi-layered cases as determined by the 35-GHz cloud radar. 
Figures 4 and 5 show scatter plots for the water and ice cases, respectively. For larger ODs, very few 
single- layer clouds occur with r e > 1 5 pm. However, some multi-layered clouds have smaller droplet 
radii. Both large and small droplet sizes occur for the smaller ODs with no clear distinction between the 
single- and multi-layered clouds. Pixels that are only partially cloud filled could be responsible for the 
larger values of r e at small ODs. In the ice c;tse, the multi-layer clouds generally correspond to < 

70 pm for x > 6 . At smaller ODs, large and small ice crystals occur for both single- and multi-layered 
systems. Histograms of the particle sizes for x > 1 0 are shown in Figures 6 and 7. As noted before, the 
largest droplets (Figure 6 ) correspond primarily to overlapped clouds, but many of the multi-layer 
systems yield values of r e that are within the observed probability distribution for single -layer clouds. 
Better discrimination may be possible for the ice clouds because the single- and multi-layered D e 
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tive droplet radius for single- and multi-layered liquid phase clouds from GOES-8. 
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Figure 5. Same as Figure 4, except for ice phase only. 
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Figure 7. Same as Figure 6, except for cloud ice crystal diameter. 
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histograms are substantially different with peaks at 50 and 80 pm, respectively. Using a threshold of D e 
= 70 pm would identify 84 percent of the multi-layered clouds having a total OD greater than 8. 
However, it would misidentify 23 percent of the single-layered clouds having relatively small ice 
particle sizes thus, skewing the ice crystal size distribution. By combining both methods, it may be 
possible to account for some of the weaknesses in the two techniques used alone. For example, the 
BTD45 data provide minimal discrimination for overlapped clouds identified as ice, but the particle size 
can be used to detect many of the overlapped ice clouds. Conversely, the particle size is not particularly 
helpful in the identification of many overlapped water clouds, but BTD45 appears to be useful in this 
case. 

Discussion and Concluding Remarks 

The results found here indicate that a multi-layered detection algorithm using BTD45 with the derived 
phase and particle sizes may provide accuracies of ~80 percent for cloud systems having optical depths 
greater than 8 to 1 0. Detection of multi-layered systems with smaller ODs or with better accuracy 
requires much additional study. The radar cloud multi- layering classification used a broad definition of 
cloud overlap. Precipitating clouds, broken clouds, and water- over- water and ice-over- ice clouds were 
included in the dataset. The amount of cloud layer separation was not specified, so that two layers only 
several hundred meters apart may have been included. By further subsetting the dataset using this 
additional information, it should be possible to refine the criteria needed to determine the presence of 
multi-level clouds. It will also be possible to define the conditions when this method is applicable. 

This paper is the first step in developing a robust method for detecting multi-layered clouds using the 
satellite imager channels available for the most routine monitoring of clouds over the ARM sites. 

Further study will be devoted to examining the impact of cloud temperature, fraction, and layer gaps on 
the thresholds that could be used for identifying multi-layered clouds. By gathering similar datasets 
over the NSA and TWP sites, it should be possible to adjust any algorithm for the particular 
environment. 
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